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ABSTRACT; 



In this paper the standard wing equation, as normally derived from 
lifting line theory, is further refined and a solution procedure more 
basic than the usual collocation technique is developed. The calculation 
method adopted avoids the necessity of performing an explicit matrix 
inversion; all equations can be solved sequentially, one at a time. On 
the other hand this technique involves the evaluation of numerous integrals 
over the span. The calculations are cumulative, and can be carried as far 
as necessary to achieve any required degree of accuracy. The analysis is 
interesting not only for purposes of practical calculation but also for 
the light it sheds on the essential mathematical structure of the basic 
aerodynamic phenomena involved. This same general method of calculation 
can also be readily adapted to the solution of other common types of 
engineering problems. 
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1. INTRODUCTION 



The standard wing equation, as usually derived on the basis of ele- 
mentary lifting line theory, is normally solved by the collocation technique. 
While this approach gives acceptable results under normal conditions, 
circumstances can arise which make a somewhat deeper analysis desirable. 
The first purpose of this paper is to show that the standard wing equation 
itself can with advantage be further developed and refined. In particular, 
it is profitable to split this equation into two parts which respectively 
govern the basic and additional lift distributions, and to refine the result- 
ing pair of equations by an integration process. The second purpose of 
the paper is to develop a calculation technique which is somewhat more 
fundamental and powerful than the usual collocation procedure. The method 
selected is based on the theory of orthogonal functions. It provides an 
explicit solution procedure which can be carried out as far as necessary 
to satisfy the initial equations to any required level of accuracy. The 
solution obtained is equivalent to that obtainable by matrix inversion, but 
the present technique achieves this result by an indirect process in which 
an explicit matrix inversion is not required at any point. The development 
of this method is of interest not only for purposes of practical calculation 
but also for the theoretical light it sheds on the mathematical structure of 
wing aerodynamics. The same general calculation technique can also be 
readily adapted to the calculation of other types of engineering problems 
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as well; the wing problem represents only one particular application of 
the method. 

All symbols are defined in Section 2 . For readers interested primarily 
in the resulting method of calculation, rather than in the detailed analysis 
itself, the calculation procedure is summarized in Section 3. The essen- 
tials of the standard wing equation are reviewed briefly in Section 4 and 
the suggested refinements are presented in Section 5. Solution by the 
usual collocation method is briefly discussed in Section 6. Orthogonal 
functions, which provide the basis for the proposed alternative method of 
solution, are analyzed in Section 7, and the resulting solution is developed 
in Section 8. 

2. DEFINITIONS OF SYMBOLS 

dC, 

a wing lift curve slope , f 

da 

a , coefficients which occur in Eq. (7.8) 
nk 

a section lift curve slope 

00 

a mean section lift curve slope as defined in Eq. (3.4) 

00 

A initially unknown constants which define the additional lift 
distribution, Eq. (3.34) 

A* constants defined by Eq. (8.11) 

AR aspect ratio 2 s/c 

b^j^ elements in orthogonality matrix, Eq. (3.12) 

B initially unknown constants which define the basic lift distribution, 

^ Eq. (3.33) 



2 



c local wing chord at spanwise station Tl, ft. 

c mean chord, S/2s, ft. 

c section lift coefficient at location Tl 
I 

C initially unknown constants which define the circulation function, 

^ Eq. (4.3) 

wing lift coefficient 

i-i 

induced drag coefficient of wing 
i 

f forcing function defined in Eqs. (3.23) or (3.24) 

forcing constants defined in Eq. (3.25) 

I planform integrals defined in Eq. (3.9) 

n 

wing integral defined in Eq (3.1 1) 

integrals defined in Eq. (7.6) 

s wing semi-span, ft. 

, 2 

S wing area , ft . 

V remote velocity, ft/sec 

00 

w downwash velocity at wing station T], ft/sec 

O' angle of attack of wing, that is, angle between zero lift direction 
of mid- span section and V^, rad 

a angle of attack of wing at which resultant wing lift vanishes, rad 

Ij 

o 

angle of twist, that is, angle between zero lift direction at station 
T| and zero lift direction at mid-span, a prescribed function of T) 
or 0 , rad. 

O' twist angle at wing tip, rad 

h 

2 , 

r circulation about wing at spanwise station Tj, ft /sec 
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6 V Induced drag factors defined in Eqs. (3.30), (3.31) and (3.32) 
BA ’ 

6 ; 

A 

6 , Kronecker delta , defined in Eq. (7.2) 
nk 

T] spanwise location expressed as fraction of semi- span 
9 angular coordinate defined in Eq. (3.1) 



"b 1 

V basic and additional lift functions defined in Eqs. (5.25) and (5.26) 

K ^ 

orthogonal lift function defined in Eq. (8.23) 
u planform function defined in Eq . (3 . 6) 
a auxiliary parameter defined in Eq. (3.5) 

T twist function defined in Eq. (3.7) 

T mean twist defined in Eq. (3.8) 

T relative angle of zero lift as defined in Eq. (3.27) 

J-i 

o 



$ 

n 

Y 

n 



Y ' 
n 







orthogonal functions defined in Eqs. (7.2) and (7.4). 
wing functions defined in Eq. (3.10) 

wing functions defined in Eq. (4.14) 

basic and additional downwash functions defined in Eqs. (5.27) 
and (5.28) 



Subscripts 

i, k, n indices which may take on the values 3,5,7 N 

A pertains to additional lift 

B pertains to basic lift 
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Matrix Symbols 



column vector containing elements 3,5,7 N 

column vector containing elements 3,5,7 n < N 

row vector containing elements 3, 5, 7, -- N 

a square matrix with rows and columns corresponding to 
indices i = 3, 5, 7, N. 

A square matrix with rows and columns corresponding to 
indices i = 3 , 5 , 7 , n<N. 

the identity matrix 

3. SUMMARY OF THE CALCULATION METHOD 

The proposed method of calculation is summarized in full in this 
section. Only such explanations are included here as are essential for 
making the method of calculation clear. All questions of derivation and 
interpretation are deferred to other sections of this paper. The formulas 
are listed in the approximate order in which they would be used in the 
calculations . 

Preliminaries 

T\ = cos 9 Change of variable (3.1) 

— s 

c = “ Average chord (3.2) 

^ s 
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AR = 



2s 



Aspect ratio 



(3.3) 



n 



a = ; a 

00 00 



(f) 



dT] Mean slope of section lift curve 



(3.4) 



CT = 



7TAR 



An auxiliary parameter 



(3.5) 



u(0) = 



T (0) = 



a 


/_c_l 


a 

00 ' 


Id 



cx 



Planform function 



Twist function 
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T = 



TudT] Mean twist 



'-7T/2 



(3.6) 



(3.7) 



(3.8) 



I = 
n 



u sinn6d0 n=3,5,7, — N Planform integrals (3.9) 



Y = 
n 



nk 



4 . ^ . / s in n 0 

~ s in n 0 + na u — ; — ;r - I 

TT 1 sin 0 n 



n = 3 , 5 , 7 , N 

7T/2 

= I 

Jo 



Wing functions 



(3.10) 



Y Y d0 
n k 



n, k = 3, 5, 7 N Wing integrals (3.11) 



Orthogonality Matrix 

The matrix [b] is calculated as summarized below. It has the structure 



6 



0 



0 



0 



[b] 



33 



53 



73 



B 



N3 



55 



75 



N5 



77 



N7 



0 

0 



I 

I 




(3.12) 



The first three rows are calculated as follows 



Row n = 3 
-2 



^33 ^33 



Row n = 5 



^53 "^33 ^53 



"2 2 
^55 ~ ^55 " ^53 



*^53 *^55 ^53 ^33 



Row n = 7 



r ^ 
^'73 

1^75 



33 



53 55 



'73 

1^75 




(3.13) 

(3.14) 

(3.15) 

(3.16) 

(3.17) 

(3.18) 

(3.19) 
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Any desired additional number of rows can be calculated by repeated 



application of the following sequence, namely 
Row n 

r 



3 O 

i n3 






I 



n5 






V= 



33 



53 



a o' 
n,n-^i 



• b 



b ^ = J 

nn nn 



I a 



0 


0 1 


'„3 


^55 


0 

n 


cn 


3 ^n-2,5 ” 


U 

n-2 , n-2 




T 

'n3 ; ^n3 






' -7 a _ 

n/ n7 







(3.20) 



(3.21) 



la oi a 
V n,n-2j n,n-2; 



fb 



r 



n3 



n5 



1= b 



n3 



n5 



nn 



V ^n,n-: 



i Q 

n,n-2 j 



^33 “ 



53 55 



^n-2, 3 ’^n-2,5 



0 

0 






1(3.22) 



n-2 , n-2 



Basic and Additional Lift 



To calculate the additional lift distribution, set 



f(6) = u (9) - “ sin 0 



Additional forcing function (3.23) 
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To calculate the basic lift distribution, set 



f (0 ) = u (0 ) [t (6 ) - T ] Basic forcing function 



(3.24) 



and change the constants A below to B . 

n n 



^/2 

F = I f(0) ijf (0) d0 
n 1 n 



J 



Basic or additional forcing (3.25) 

constants 



n = 3, 5, 7, — N 



Solution Constants 

^ A '' 

^3 



< *7 






5 . 



>= [b] [b] r 

I 

—1 






VF 



N 



J 



(3.26) 



Wing Characteristics 

N 



°Lo 

'"t. 



= T, 



Y' nb 1 - T 

/ n n 



nB I 
n r 

n = 3 , 5 , 7 , 



Zero lift angle 



(3.27) 



N 



a • 

00 



a / 



i = 1 +a 



( ' ) 



Wing lift curve slope 



n = 3 , 5 , 7 



(3.28) 



D. 



nm I CO t^ 






B 



(3.29) 



Induced Drag 
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where 



N 

n=3 , 5 , 7 , — 
N 

= V nB A 
BA ) n n 

n=3 , 5 , 7 , — 

N 

2 

b = r\h 
A ) n 

n=3 , 5 , 7 , — 



i 

I 

} 

i 



\ 



I 

} 



I 

t 



1 



/ 



Induced 
drag 
factors . 



Lift Distribution 



N 



c 

I B / _c_ 

a" a — ' — 

“ t^ 1 C / 



\ B sin n 9 

L n 



n = 3 , 5 , 7 , — 



Basic 




4 

7T 



N 

sin 0 + V A sin T1 G 

z ^ 

n=3 ,5,7, — 



Additional 



(3.30) 



(3.31) 



(3.32) 



(3.33) 



(3.34) 



4. REVIEW OF THE STANDARD WING EQUATION 

By utilizing the standard lifting line idealization of a wing, and em- 
ploying the Biot Savart law, we can deduce [1 , 2] * that the downwash w 
at any point T] along the span is related to the distribution of circulation 
r along the span as follows: 



* Numbers indicate references listed at end of paper. 
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w 



V 

00 



1 

TT 



+ 1 



J 

-1 



dll 



4V s 



dTl' 



(T)' - T\) 



(4.1) 



It is customary to introduce the change of variable 



T) = cos 9 
T|' = cos 9’ 



(4.2) 



and to describe the distribution of the circulation Fby a series of the form 



^ C„slnne’ (4.3) 

\ 00 * 

n = 1 , 3, 5, — 

The present discussion is restricted to lift distributions which are 
symmetrical about the mid-span and hence only terms of odd index are 
retained in Eq. (4,3). 

Upon introducing Eq. (4.3) into Eq. (4.1) and making use of the 
integral formula 



1 

IT 



77 

r 



cos n 9' d 9‘ 
cos 9' - cos 9 



sin n 9 
sin 9 



(4.4) 



the following solution is obtained for the downwash 



I w _ V n n sin n 9 

1 V , 2 n sin 9 

00 

n =1 , 3 , 5 , 



(4.5) 
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Once the above integration has been performed, there is no further 
need to distinguish between the variables T] and T]' or 0 and 0' and the 
prime marks may be deleted from Eqs. (4.2) and (4.3). 

The lift exerted by an element of the wing lying between Tland (T| + dT]) 
may be expressed firstly in terms of the section lift coefficient c^^ and 
secondly in terms of the circulation F, in accordance with the Kutta 
Joukowski law. Equating these two alternative expressions for the local 
lift gives 



c 



c; 



= 4 m 



-i-1 

4V s / 

OO ' 



(4.6) 



where 



c = 7“ “ average chord 

2s 



(4.7) 



2s 

■AR = — = aspect ratio 
c 



(4.8) 



The local lift coefficient can be expressed in terms of the local 

section lift curve slope a^ and the local net angle of attack, taking into 

( ■ w \ 

— jin reducing the effec- 

00 ' 

tive angle of attack. Hence we obtain 



W \ 

(a + - — ) (4.9) 

00 

where a = wing angle of attack = angle of attack of the zero lift 
direction of the airfoil section at midspan with respect 

to the remote velocity V . 

00 



12 



and a (6) = angle of twist = angle of zero lift direction at Tlwith respect 

to zero lift direction at midspan, a prescribed function of B. 

Upon combining Eqs. (4.3), (4.5), (4.6) and (4.9) and rearranging 
we can finally obtain the fundamental wing equation in the following form, 
namely , 

y C ilf ' (9) = ^ u (a + a ) (4.10) 

n n t 

n = 1,3,5, 

where 

'1 

a 1 a dT| = mean section lift curve slope (4.11) 

“ = I “lei 

I 

0 



a 



a 



7TAR 



an auxiliary parameter 



(4.12) 



1 ^( 9 ) 




a prescribed planform function 



(4.13) 



r (9) = - 



sin n P + au 



n sin n 6 
sin fi 



(4.14) 



Apart from minor changes in notation, Eq. (4.10) is the fundamental 
wing equation as found in most elementary texts on aerodynamics [1,2]. 
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5 . ALTERNATIVE EQUATIONS FOR BASIC AND ADDITIONAL LIFT 



By substituting Eq. (4.3) into (4.6) multiplying through by dT) = 

- sin 0 d 0, integrating over the semi-span and taking advantage of the 
orthogonality relation 



This verifies the well known fact that the overall lift of the wing is 
fixed solely by the first term of the series in Eq. (4.3). The other terms 
affect the form of the lift distribution but not the resultant net lift. 

It is useful to express the twist distribution of the wing in dimension- 
less form as follows 



The basic wing relation, Eq. (4.14) , can be further developed by in- 
troducing the substitution 



tt/2 



sin n 9 sin k 9 d 9 = 6 



nk 



+ 1 if n = k 



0 



0 if n k 



(5.1) 



we obtain a solution for the overall lift of the wing in the form 
.1 



C 



L 




1 



(5.2) 




a prescribed twist function 



(5.3) 



where cy^ = angle of twist at wing tip 



(5.4) 
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1 

irm 



~ 



(5.5) 



C 

n 



a a B + A 
”t^ n L n 



n = 1 , 3 , 5 , 



The constants B now define the basic lift distribution of the wing, 
n 

This is the distribution associated with the wing twist. Specifically it 

is the distribution which occurs when the wing is at the angle of attack 

O' = O' which gives = 0. The constants An define the additional 
o 

lift distribution that occurs as a result of changes in angle of attack. 

Hence it is proportional to the angle {a - a. ) or to itself. 

Lo L 

Consider Eg. (5.5) for the specific case n = 1. Substitute for 
from Eg. (5.2) . The resulting expression must be satisfied for all values 
of C. . In particular by setting = 0 we find that 




A^ = 1 (5.6) 

Substituting the expressions (5.5) into the basic wing relation Eg. 
(4.10), and making use of the special result (5.6), we obtain 






= cr u ( 0 -+ a) (5.7) 



n = 3 , 5 , 7 , 



n= 3,5,7, 



This must be satisfied for all operating conditions. In particular for 
O' = and = 0, it reduces to 
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(5.8) 



a»y s'*" (6) = 0‘u (» + aj 

t.Z- n n t 



n = 3 , 5 , 7 , 



Then subtracting Eq. (5.8) from (5.7) we obtain 



Sin6+y 



A Y' (e) 
n n 



= a u ( O' - ) 



(5.9) 



n = 3 , 5 , 7 , 



Eq. (5.8) can be further reduced by dividing through by Oa and 



setting 



O' / LO 
1 ' 



Eq. (5.9) can be further reduced by substituting 



(5.10) 






(5.11) 



( dC 

L) = wing lift curve slope 

da 



(5.12) 



In this way we obtain for the basic and additional lift the two funda- 
mental equations. 



Y B f (0) = u (t + t) 
; n n Lp 

n = 3,5, 7 , _ 

sin 0 +y A y' (e) = ( — I 
Ann a / 



n = 3,5,7, 



(5.13) 



(5.14) 
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The complete analysis of a twisted wing can be obtained either by 



solving the standard wing relation, Eq. (4.10) for two different angles of 
attack, or alternatively, by solving the above pair of relations, Eqs. (5.13) 
and (5.14), for basic and additional lift. The total number of unknowns to 
be found is exactly the same by both methods; only the form is different. 

It can be said, however, that Eqs. (5. 13) and (5.14) are somewhat more 
fundamental in form in the sense that, for a wing of prescribed geometry, 
the unknowns B , t„, A and! are true invariants, whereas each one 

" « " \rl 

of the "constants" C in Eq. (4.10) is actually a function of angle of 

n 

attack . 

Moreover, Eqs. (5.13) and (5.14) lend themselves to further develop- 
ment in a useful and interesting way. In this connection, the first step 
is to multiply these equations through by dTj = -sin 0d9 and integrate 
over the semispan. 

This operation involves the following integrals. Firstly, 



1 

udTl = 

0 




(5.15) 



a result which follows directly from the definition of a as given in 

00 

Eq. (4.11). 



Secondly , 
1 



tr/2 



4 

7T 



sin n© dt? = ~ 
7T 



sin n6 sinSdG - 1 for n = 1 

= 0 for n ^ 1 (5.16) 



J 

0 
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Thirdly, we define the average value t of the twist function by the 



integral 



ir/2 

f 



= \ TudTl = t TV sin 9de 



J 



(5.17) 



Finally, we define the series of integrals 

1 ir/2 



sin n 0 ^ 

I = \ V — ; — - — dll = \ V sin n 9d 0 

n sin 0 



J 



(5.18) 



n = 1,3,5, 



We also note from this last definition that specifically for n = 1 



1 

f 

1 udTl = 1 



(5.19) 



Upon integrating Eqs. (5.13) and (5.14) over the semi-span, utilizing 
the foregoing integrals, then simplifying, we obtain the two significant 
results below, namely. 



T 







nB I 
n n 



T 



n = 3,5,7, 



(5.20) 
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(5.21) 



— = 1 + a 

a 



1 +y nA I 
I L n n 



n = 3,5,7, 



These results are useful for two reasons. First of all, they can be 



used to eliminate the unknowns t and[_^\ from Eqs. (5.13) and (5.14) 

^0 «a 1 

thus reducing the number of simultaneous unknowns which must be found 
from these equations. Secondly, they provide a subsequent solution for 
these two unknowns, and one that is especially accurate because it is 
based on an integration over the entire wing. 



After Eqs. (5.20) and (5.21) have been used to eliminate andi 

■"o 



( \ 



a I 



from Eqs. (5.13) and (5.14), and after the results have been simplified, 
they become finally 



y 



Kt (0) = u (t - t) 
n n 



(5.22) 



n = 3,5,7, 



A ilr (0) = (u - - sin ) 
n n IT 



(5.23) 



n 3,5,7, — — — — 
where 

- I^j (5.24) 

n = 3,5,7, 

Eqs. (5.22) and (5.23) represent an alternative form of the funda- 
mental wing equations. They can advantageously replace the standard 



, 4 . . ^in n9 

I’ (0) T sin n 0 + nau| — — — 
n sin 0 
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wing equation (4.10) as usually employed in the lifting line analysis. 
The series on the left side of Eqs. (4.10) and (5.22) and (5.23) are, of 



course, truncated at some value n = N, depending on the accuracy 



Eq. (4.10). Note that in this case, Eqs. (5.22) and (5.23) will contain 



over the latter solutions can be expected to be somewhat more accurate 
since they include the averaging effects of integration over the semi-span. 

Note that Eqs. (5.22) and (5.23) are of identical format on the left. 
Only the forcing functions on the right are different. Eq. (5.22) shows 
clearly that the constants B depend only on the dimensionless planform 



and twist functions and, in particular, that an untwisted wing gives 

= 0. Likewise, Eq. (5.23) shows that the constants depend only 

on the deviation of the planform function u from an ellipse of equal area 

and, in particular, that a wing having an elliptical planform function 

gives A = 0 . 
n 

It should also be noted that while Eqs. (5.22) and (5.23) are more 
effective than Eq. (4.10) in the sense noted above, they entail the addi- 
tional labor of evaluating the integrals t and 

After Eqs. (5.22) and (5.23) have been solved for the constants B^ 



required. Often only the four terms in , C„, C and C are retained in 

-1- O O / 



only three unknowns apiece, namely, B. , B , B and A , A . A . More- 

o O / o o / 



n 




T and the slope of the wing 



lift curve 




can then be found accurately from Eqs. (5.20) and (5.21) , 



respectively. 
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The resulting basic and additional lift and downwash distributions can 
be found by combining Eqs . (4.3), (4.5), (4.6) and (5.5). We thereby 
obtain the following solutions. 



ic 



XB \ 



c \ 



- ii- 

a a / c 

“h" 



^ y B sin n Q = \„ (p) 
ir Ls n B 

n = 3,5,7 , 



(5.25) 




(5.26) 



(5.27) 








n sin n 6 
n sin 0 

n = 3,5, 7 , 




(&) 



(5.28) 



Once the constants A and B are known, they fix the four functions 

n n 

above which define the fundamental aerodynamics of the particular wing 
under consideration. 

The one remaining question of basic significance which falls within 
the purview of lifting line theory concerns the distribution and overall 
value of the induced drag. The overall coefficient of induced drag is 
defined by the integral 

1 J- 



C 



Di 




(5.29) 
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Upon substituting for the quantities in Eq. (5.29) in terms of the 
four functions defined in Eqs. (5.25) through (5.28), then rearranging, we 



obtain 



n/2 



,2 



a a 
1 ^ t 

^ 

i AR 



Xg (9) (9) sin 0d9 



7T/2 



+ 



a O' C. 
” t^ L 

7T AR 



ff /2 



i«A (0) + (0) a'g (0)] sin 0d0 



J 

0 



+ L 



7TAR 






sin Gd 0 



(5.30) 



The integrals on the right side of Eq. (5.30) show that there are three 
distinct components of induced drag and the respective integrands define 
the corresponding distributions. However, instead of a direct evaluation 
of these three separate integrals, it is simpler to proceed as follows. 

First we evaluate the overall integral in Eq. (5.29) using Eqs. (4.3), (4.5) 
and (4.6). The result is simply 

= -rr m ) nC (5.31) 

D. / n 

1 

n = 1,3,5 , 
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Next we substitute Eq. (5.5) into (5.31), expand and regroup terms. In 
this way we obtain finally 



IT m 






where 6 



B 



BA = 



n = 3,5,7,. 

y nB A 
/ n n 

n = 3,5,7, 

= y nA^ 

/ n 

n = 3,5,7, 



(5.33) 



(5.34) 



(5.35) 



We also note in passing that for the special case of an aerodynamically 
untwisted wing of elliptical planform, Eq. (5.26), (5.28), (5.21), and 
(5.32) reduce, respectively, to the following familiar and elementary 
forms, namely. 



— I = ~ sin 0 = elliptical lift distribution 



(5.36) 



w 

V 



ITfR 



= uniform downwash along span 



(5.37) 



a = 



1 + ^°° 
1TAR 



— = slope of wing lift curve (a maximum) 



(5.38) 
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induced drag (a minimum) 



(5.39) 




1 



These last four results illustrate the fact that the shift from the 
standard form (4.10) of the wing equation to the revised forms (5.22) and 
(5.23) does not affect the essential character of the solution obtained. 

6. THE COLLOCATION METHOD 

The usual method of obtaining an approximate solution of the standard 

wing relation, Eq. (4.10) is by collocation, that is, by satisfying this 

equation exactly only at a limited number of discrete points along the 

semi-span. The number of points must of course equal the number of 

unknown constants Cj^, , C^, to be retained in the approximate 

solution. The location of the collocation points is somewhat arbitrary, 

but it is customary to space them at equal intervals of 6 along the semi- 

n 

span, to include the midspan point at 9~”^/ ^rid to omit the wingtip point 
0=0. The solution of the resulting set of simultaneous equations entails 
a matrix inversion, but is otherwise quite straightforward. If a computer 
be utilized, of course, even the matrix inversion poses no appreciable 
computational burden. 

In practice the above calculations even if performed only on a desk 
calculator, usually work out in a reasonably satisfactory manner. This 
is partly owing to the fact that any well designed wing of moderate to high 
aspect ratio has a lift distribution that does not deviate strongly from the 
optimum elliptical distribution. Consequently the successive constants 
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C , C„, C , typically diminish rapidly in magnitude, and only a few 

1 o 0 

terms, often just four, suffice to produce results of an accuracy adequate 
for many engineering purposes. Moreover, a quest for extreme numerical 
accuracy would not be realistic because the basic lifting line idealization 
on which the wing equation is based is itself an approximation. 

Nevertheless, there are situations in which the mathematical circum- 
stances are not so favorable, for example, in cases of strong discontinuities 
in the lift distribution owing, say, to the action of wing flaps over part of 
the span. Moreover, apart from practical questions of computation, col- 
location has a number of features which are not optimal from the theoretical 
viewpoint. For example, the choice of the actual collocation points in- 
troduces an element of arbitrariness into the solution. It would be in 
some respects preferable to have a method whose accuracy depends in a 
known way only on the number of terms retained in the solution, and not 
on such contingencies as the location of collocation points. Furthermore, 
the collocation method can under certain adverse circumstances lead to a 
set of equations which are numerically ill conditioned. Moreover, the 
purely theoretical conception of how the collocation solution approaches 
its limit as the number of terms retained in the solution hypothetically 
increases without bound tends to be somewhat obscure; it entails the 
conceptual difficulty of visualizing the inversion of a matrix in the limit 
as the number of rows and columns increases without bound. 
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7. ORTHOGONAL FUNCTIONS 



There is an alternative to the collocation method that deals with the 

above difficulties, both conceptual and computational, in a very clear and 

decisive way. It can be developed through the introduction of a certain 

family of orthogonal functions f (9), n = 3,5,7, that are related in 

n 

a simple way to the basic functions (9) that occur in the wing equation. 

The mathematical structure of these auxiliary functions is developed in 

this section. It turns out that the functions § , while essential to the 

n 

theoretical development, can be eliminated from the final calculations; 

only certain constants that occur in the definition of these functions carry 

over into the resulting calculation procedure. 

The auxiliary functions f (9) are related to the original functions 

n 

ijf (9) by equations of the following form: 
n 



' ^3 +3 



f(9) = b 'll + b 'll 
5 ^ ’ 53 ^3 55 ^5 



f (9) = b 'll + b 'll + b 'll 
7 ^ ‘ 73 ^3 75 '^5 77 ^7 



(7.1a) 



'n(') = ^3 ^3 ^5 ^5 



+ b 'll 

nn n 



(7.1b) 



The constant coefficients b , in these equations are chosen so that 

nk 



the functions f (9) form an ortho-normal set, that is, so that 
n 
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7T/2 



$ $ d9 = 6 =0 fork = 3,5,7, (n-2) (7.2a) 

n k nk 



= +1 for k = n 



(7.2b) 



It is convenient to develop much of the theory in matrix format. Ac- 
cordingly, we introduce the following abbreviated matrix notation. 



( 



I 



n 



{ ^3 (e) 

^5 



< ^7 (e) ; 

f 



. (e) 

V. n 



n 



tb]n = 



33 



53 



73 



n3 



55 



75 



77 



' $3(9) ■ 

$5 (e) •. 

" $7(6) ^ 



' $ (e) I 

\ n y 



0 

0 

0 



b b 

n5 n7 



nn 



Consequently, Eqs. (7.1) may be written simply as 

(; 



= [bi„ ♦ 



• n 



Jn 



(7.3) 



(7.4) 



where we call [b] the orthogonality matrix. 

In connection with the ensuing analysis, it will also be convenient 
to define the following auxiliary integrals, namely. 
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(7.5) 



Jr/2 



^nk ^kn 



t (e) t (0) de 

n k 



and 



K 



hk 



n/2 

r 



ilr (0) $, (0) d0 
n k 



k = 3,5,7,— -(n-2) 



(7.6) 



The integrals will ultimately be eliminated from the analysis. 

The integrals however, can be assembled to form the following matrix. 







hs 


^37 


hn 




'53 


^55 


' — 1 
cn 

1 

1 


Jsn 


T 










n n 


CO 
1 \ 


^75 


J 77 


^7n 




i ^n3 




N _ — 

n7 


f„n 



(7.7) 



It can be seen that the matrix [b] „ fixes the functions 

n- Z o 5 / 

§ Now suppose that this matrix is known and that we wish to add 

n- 2 

t h 

to it the initially unknown n row of coefficients which will define the 
next function For this purpose, it is useful to write temporarily in 
the alternative form 



n-2 



$ = b 

n nn 



Y a $ + t|r 

2 nk k n 

k = 3,5,7, 



(7.8) 
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In matrix notation, this may be written in either of the two equivalent 



forms below, namely. 




Also, from (7.4) 




Now premultiplying Eq. (7.10) through 




(7.9) 



(7.10) 



(7.11) 



(7.12) 



(7.13) 



We now integrate this over the semi-span, that is, from 9 = 0 to 

_ E 

2' then apply the ortho-normal condition (7.2) and the definitiions (7o5) 
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and (7.6). The following simplifications occur, namely, 



f de = 0 
n 




(7.14) 



de 



[I] 



n-2 



n-2 



1 


0 


0 — 


ol 

i 


0 


1 


0 — 


0 i 

i 


0 


0 


1 — 


0 


0 


0 


0 — 


1 



(7.15) 



where [I] is the identity matrix, with rows and columns corresponding 
n ^ 

to indices 3, 5, 7, (n-2). 



Also 

7t/2 



r 

I 



$ > A d9 
• n 

n-2 



i 1 

tK \ 
1 n 



K 

{ n5 



n-2 



[k„3 

] 

< K 

t 

1 



h7 1 

- 1 

^n,n-2j 



(7.16) 



With the aid of relations (7.14) (7.15) and (7.16), the integrated 
form of Eq. (7.13) can be rearranged to yield the following significant result. 



[I] 



n-2 



f ^ 


i ^ 


i 1 


a ; 


I- = ia } 


II 


1 ^ 


1 ^ 




L J 


L_o 1 


i J 



(7.17) 



n-2 



This defines the solution for the constants a , such as to satisfy the 

nk 

orthogonality conditions (7.2a) . However, a further reduction of the right 
side of Eq. (7.17) is possible and useful. Note that, from the definition 
(7.6) and from (7.4), 
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7T/2 



K 



n 






n-2 



I 

j , 

n-2 



de = 



= [b] A 

n-2 



A 

J 



n 



= [b] 



n-2 



n-2 



f ] ' 

J ""’n-2l ♦] 



de 

^ - n-2 



> 

J. 



; 



n3 

n5 

n7 



V. ^n,n-2y 



(7.18) 



Combining (7.17) and (7.18) gives the orthogonality relation in the 
desired final form, namely, 

( 



( ■ "='n-2 1 ^nl 



(7.19) 






n-2 



n-2 



To apply the normalizing conditions (7.2b), we first post-multiply 

-2 

(7.9) by (7.10), then multiply through by b . This gives 

nn 





1 1 




-2 2 i 


1 1 i 




II 

1 

XI 


a > y 


nn n 


l"j. 


l j 



f r [ 



T f ' 






< a 



n 



- T 



+ ^ 
n t 



n 



+ I 



n 



C ^n- 



n-2 



n-2 



f 1 T ' 

i J * 

+ •( a y ‘'.$.■'111 



"J ^ 



J 



n 



n-2 n-2 



(7.20) 



We next integrate this over the semi- span noting that, from (7.2b), 
77/2 



e de = 1 

n 



(7.21) 
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and from (7.5) 



tt/2 



i 

I 

I 

J 



0 



2 

t d9 
n 



nn 



(7.22) 



Also, the first term on the right side of (7.20) may be simplified after 
integration in the manner of (7.15). The next two terms, which incidentally 
are of equal magnitude, can be simplified according to (7.16). In this 
way we ultimately obtain from the integrated form of (7.20) the result 





+ i K 



+ J 



nn 



n-2 



n-2 



(7.23) 



This reduces readily to 



f 



b-2 

nn 



i 

V ^ - V 



v_ 



n 



i a 



n 



J 

n-2 



J 



-T f 



+ 2 ' a 
' n 



K ! 
n 



+ J 



n-2 



I i u 

n-2 n-2 



nn 



(7.24) 



We can eliminate / K 

r 

tion in the form given in (' 



from this by using the orthogonality rela- 
The final result therefore becomes simply 
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b 



-2 



nn 




(7.25) 



At this point the three relations (7.10), (7.19) and (7.25) together 
define the function in such a way that it satisfies the normalizing con- 
dition (7.2b) and is orthogonal to all the previous functions in the set, 

namely, to f_, $ However, it still remains to determine 

o 5 / n-2 

the coefficients b ., b _, b _ b . which are needed to describe 

rio n5 n / n / n ^ 

$ in the standard format of Eq. (7.4) . In this connection we note first 
n 

that 



$ 

n 




n-2 n-2 



1 

♦ 



n 



(7.26) 



Now a comparison of (7.26) with (7.9) shows that 




(7.27) 
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Moreover, by analogy with Eq. (7.4) we may write 
^ r 



‘ $ ' 

'' n-2 



[b] 

^ n-2 ^ ^ j n-2 



(7.28) 



Substitution of (7.28) into (7.27) discloses finally that 
T T 



■^b'^ = b J a ^ 

n i nn ] n ' 

i •' ? 

^ n-2 n-2 



[b] 



n-2 



(7.29) 



This is the required result which, along with (7.19) and (7.25), 

th 

enables the n row of matrix [b] to be calculated when all the previous 

n 

rows are known. 

The principal results of this section can now be summarized. It is 

supposed that all quantities over the range of the index i = 3,5,7, 

(n-2) are known from previous calculation. The following calculations 

then extend these results to the additional index i = n. 

Firstly, we must evaluate the additional integrals 

n/2 




r 

(6) ♦ (e)d(i 

n k 



k = 3,5,7, n 



(7.29) 



J 



0 



For any given wing planform, the functions as defined by Eqs. 

(5.24) are known. Hence the constants J , are computable by direct 

nk 

numerical integration. 
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Next we calculate, in the order given, the quantities 



. a 

n 




n-2 



T 



b~^ = J - • a 
nn nn n 




(7.30) 



(7.31) 



■nT 



' b 



n 



n-2 

( 

b , a ’ 
nn / n / 



n-2 




(7.32) 



n-2 \ n-2 

Eqs. (7.31) and (7.32) fix row n of the orthogonality matrix [b] . 

n 

Consequently, the function $ is now defined in the standard format, namely, 

n 

T 

$ =' b < Y (7.33) 

n n 



n n 

Moreover, the function f so defined satisfies the ortho-normal rela- 

n 

tion (7.2) with respect to all previous functions in the set. 

The above calculation sequence, Eq. (7.29) through (7.32), may then 
be repeated for the next successive row of the matrix [b] . In this way, 
the fundamental matrix [b] can be expanded to any desired number of rows 
and columns . 

We shall next show how, once the orthogonality matrix has been es- 
tablished to any desired size, all remaining unknowns of the problem can 
be calculated to the corresponding degree of accuracy in a routine way. 
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8. SOLUTION BY MEANS OF THE ORTHOGONALITY MATRIX 



Both of the revised wing relations, Eqs. (5.22) and (5.23), are solved 
by the same general method and hence for the purposes of this section, it 
is preferable to adopt a generalized notation which embraces both equations. 
Also, it is helpful to retain the matrix notation. Hence we rewrite Eqs. 

(5 .,22) and/or (5.23) simply as 



< ♦ f ' A I 



( J 

where 



t3 (9) 

If .(9) 



J ^ > = J ^ (9) " 

I i » / 



I J 






f (9) 



and 



A ^ 



^ ' 



N 



( 8 . 1 ) 



( 8 . 2 ) 



For calculating the additional lift distribution we set 



f(9) = f^(e) 



o - “I sin 0 
7T 



(8.3) 



For calculating the basic lift distribution we set 



f(9) = fg(9) = V (t- t) 



(8.4) 



In this latter case we also change j a| to | Bjfor consistency with the 
notation established earlier. 



It is also convenient to define the set of integral constants 

n/2 



F = 
n 



f(9) t (9) d0 
n 



n = 3,5,7, — N 



(8.5) 
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and, of course, the corresponding vector 



F ^ 

■ / 

1 



< F 



= / 

> 



( 8 . 6 ) 



\J 



Also since 

r 1 • 

i 5"= [b]^ r 

it follows that 






(8.7) 



( 8 . 8 ) 



and that 



T T 

1*1 ={*) 



[b] 



and that 






-T 



(8.9) 



[b] (8.10) 

This last result suggests the usefulness of defining the auxiliary 



vector 



-T 



A* |- = [b] •( A 



J 



( 8 . 11 ) 



We now find that the basic relation to be solved. Eg. (8.1), can be 



rewritten in the alternative form 

.T 



$ 



5i-) 



f(9) 



( 8 . 12 ) 



The great advantage of (8.12) over the original form (8.1) lies in the 
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orthogonality properties of the functions For if we now multiply (8.12) 

through by any one of the functions $ (6) , then integrate over the semi- 

n 

span, we find that, because of the orthogonality properties, all integrals 

but the one in n vanish on the left with the simple result that 

7r/2 



A* = 
n 



) f(e) $ (P)d 9 



n 



(8. 13) 



n = 3,5,7, (n-2) 



Upon substituting into (8.13) the relation 



,T . 
1 ! 

5 



f = ^ b 



n 



n 



1 



I 



n 



n 



(8.14) 






and utilizing the definitions (8.5) , this result can be rewritten in a form 
more convenient for calculation, namely. 




(8.15) 



This can readily be generalized over the entire vector 
process amounts to pre-multiplying Eg. (8.12) through by (8.7), then in- 
tegrating over the semi-span and applying the ortho-normal conditions in 



the usual way. We thereby obtain 

i' ; f \ 

I A* = [b] \ F ^ 

] I 



(8.16) 
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A significant property of the constants A* as revealed by Eqs. (8.15) 

n 

or (8.16) is that they are true invariants of the solution in the following 

sense. The value of A* depends on F. , F_,F_, F . But owing to the 

n o 5 7 n 

zeros in [b] , it is independent of all the higher order terms F F 

n + 2 n + 4 , 

and so on. Consequently, adding on additional terms to the solution in no 

wise affects the values of the constants A* previously calculated. 

, n 

[ *' 

Nevertheless, while ^ A ’ contains the true invariants of the solution, 

✓ 

it is more convenient in some ways to continue to express the solution in 

terms of the original unknown constants, as contained in the vector A j-. 

( ■ 

To solve for i A^ we need merely to premultiply Eq. (8.11) through by 
T 

[b] . The result is simply 






= [b]' 



A* 



) 



(8.17) 



In many cases it may be preferred to bypass the vector -i A*f altogether 

i J 

in the solution. To show how this can be done, we simply substitute for 
I A*J from (8.16) into (8.17). The result becomes, finally. 



|a;= [b]^[b]|Fj 



(8.18) 



Thus, once the orthogonality matrix [b] is known, Eqs. (8.5) and (8.18) 

summarize the solution for the initially unknown constants A . When these 

n 

are fixed, all remaining features of the solution can be established from 
the equations developed in the earlier sections. 

The solution represented by Eq. (8.18), although it entails no explicit 
matrix inversion at any point, is in fact equivalent to a matrix inversion. 
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To show this, we premultiply the basic relation, Eq. (8.1), through by 
<'1'^ then integrate over the semi-span and apply the definitions of 
and as given in Eqs. (7.5) and (8.5) respectively. In this way 
we obtain 

[J] " A f = |F j (8.19) 

The solution of this equation is 

i A .' = rj]' F f (8.20) 

( ; I > 

Comparison of Eq. (8.18) with (8.20) discloses that 

[b] ^ [b] = [J]"^ (8.21) 

which demonstrates the sense in which calculation of matrix [b] amounts 
to the inversion of matrix [J] . 

Bearing in mind that the elements of [b] are calculated in the first place 

solely from the elements of [j] , it is clear that what we have here in effect 

is a particular method of inverting matrix [J] . But since [J] could be any 

symmetrical matrix, the method amounts to an inversion procedure applicable 

to any symmetrical matrix. While the method happened to be developed 

in this instance in terms of the orthogonal functions # , these functions 

n 

do not actually occur in the resulting calculation formulas. Hence the 
final algorithm m.ay be applied to invert any symmetrical matrix quite ir- 
respective of its origins. In this particular instance, the elements of [J] 
happen to be certain specific integrals; the same inversion procedure 



40 



should also apply to symmetrical matrices whose elements are defined in 



some other way. Hence the technique is clearly applicable not just to the 
wing problem, but to a wide range of other applications as well. 

One advantage of solving by means of the matrix [b] is that the cal- 
culations are strictly progressive and cumulative. The process of adding 
on each additional row to [b] provides additional finer detail about the 
solution without disrupting the details calculated from earlier rows. This 
becomes somewhat clearer if the solution be expressed in terms of the 

constants A* rather than in terms of the original constants A . The in- 
n n 

variance properties of the A* have been noted earlier. In contract with 



this, the constants A fluctuate each time a higher order term is added 



to the solution, and approach their true ultimate values only in the limit 
as the number of terms increases without bound. This analysis therefore 
sheds much light on the nature of the convergence process that, in prin- 
ciple, defines the exact mathematical solution of the original basic 
equation. 

The solution may therefore be viewed as the result of superimposing 
an indefinite number of modes each of which is orthogonal to all of the 
preceding modes. To illustrate this concept, consider the solution for 
the additional lift distribution, Eq. (5.26), rewritten in matrix notation, 
that is , 



n 



n 



T 



C 



c 




L 



( 



(8.22a) 
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(8.22b) 



where 




(8.22c) I 






: 




f 



sin 3 9 



4 

V 




sin 5ft I 

) 

sin 70 ! 



(8.23) 



, defined in Eq. 

(8.23) represent the successive invariant orthogonal modes, a linear 
combination of which fixes the additional lift distribution as shown in 
Eq. (8.22c) . The addition of any mode n to the solution does not affect 
the amplitudes of the modes 3, 5, 7, — (n-2) previously obtained. This 
example illustrates the theoretical concept, at any rate. Of course, for 
practical calculation, the ordinary form shown in Eq. (8.22a) would 
probably be preferred as involving somewhat less overall computation. 



( ^ (9) ] I Sin N e; 

Thus the functions 
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9. 
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on the essential mathematical structure of the basic aerodynamic phenomena 
involved. This same general method of calculation can also be readily 
adapted to the solution of other common types of engineering problems. 
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